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SECTION 1 


INTRODUCTION 

Future space missions such as the Earth Science Geostationary Platform (ESGP) 
will require highly accurate antennas with apertures that cannot be launched 
fully formed. The operational orbits are often inaccessible to manned flight 
and will involve expendible launch vehicles such as the Delta or Titan. There 
is therefore a need for completely deployable antenna reflectors of large size 
capable of efficiently handling millimeter wave electromagnetic radiation. 

The parameters for the type of mission considered herein are illustrated by 
the heavy shaded horizontal bars in Figure 1. This logarithmic plot of fre- 
quency versus aperture diameter shows the regions of interest for a large 
variety of space antenna applications, ranging from a 1500-meter-diameter 
radio telescope for low frequencies (less than 10 MHz) to a 20-meter-dlameter 
infrared telescope. For the ESGP, a major application Is the microwave radl- 
ometry at high frequencies (up to 220 GHz) for atmospheric sounding. The 
heavy lines in Figure 1 occur at peaks and windows of the absorption spectra 
and are useful for the determination of atmospheric temperature, clouds, water 
vapor and precipitation; the width of the lines denotes the bandwidth of 
interest. The aperture diameters start at 4 meters, the size which can be 
launched without folding, and range up to the size yielding a resolution at 
the Earth's surface of about 6 km. 

In the figure, only those frequency bands above 30 GHz are shown. These 
higher frequencies require a solid reflector surface, perhaps segmented or 
inflated. On the other hand, the lower frequencies can be reflected effi- 
ciently by expandable mesh surfaces. 

Almost all existing large antenna reflectors for space employ a mesh-type 
reflecting surface. Examples are shown and discussed In Reference 1, which 
deals with the various structural concepts for mesh antennas. Fortunately, 
those concepts are appropriate for creating the very large apertures required 
at the lower frequencies for good resolution. 

The emphasis of this paper Is on the structural concepts and technologies that 
are appropriate to fully automated deployment of dish-type antennas with solid 
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Figure 1 . Large space antenna requirements. 
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reflector surfaces. First the structural requirements are discussed. 
Existing concepts for fully deployable antennas are then described and 
assessed relative to the requirements. Finally, several analyses are pre- 
sented that evaluate the effects of beam steering and segmented reflector 
design on the accuracy of the antenna. 
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SECTION 2 


STRUCTURAL REQUIREMENTS 

A probable configuration for the high frequency radiometer antenna is shown in 
Figure 2. It consists of a primary reflector dish, a subreflector, and more 
than one feed system. For a structural point of view, each reflector consists 
of a reflecting surface and a structure to hold the reflecting surface in 
shape and position. In some cases, the two functions can be combined, but it 
is helpful to consider them separately. 

2.1 REFLECTING SURFACE 

Passive microwave radiometers must have very high efficiency because of the 
feebleness of the received signal. Thus, the reflecting surface must cause 
minimum loss. This requires a surface of high conductivity. The surface can 
be very thin electrically because the skin depth of the surface currents Is 
very small (much less than one micrometer). If the surface Is a grid, low 
loss requires that the grid spacing be a small fraction of the wavelength A, 
say A/50. Similarly, the surface must be smooth, with roughness less than 
A/50 for undulations having a spatial period of a half wavelength or more. 
Thus, the compliant knitted mesh that readily stows Into a small package is 
not suitable for frequencies greater than about 30 GHz. Breaks or gaps In the 
reflecting surface are acceptable If they are many wavelengths apart and If 
the large ones do not form a regular pattern. 

2.2 SUPPORTING STRUCTURE 

The supporting structure must be made sufficiently accurate, stiff, and dimen- 
sionally stable In order to meet the stringent requirements for diffraction- 
limited antenna performance. Not only must the antenna be efficient, but also 
must It exhibit small side lobes. Analysis (see Reference 2) shows that 
large-correlatlon-dlstance surface errors with an rms of A/50 can raise the 
near-in side lobes by as much as 20 dB down from the main lobe. In addition, 
any distribution of surface normal errors with an rms of A/50 will reduce the 
main-lobe efficiency by six percent. It appears, therefore, that a demanding 
mission such as microwave radlometry requires a smaller rms error, probably 
A/100. 
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Figure 2. Example of multiband high-frequency radiometer antenna. 
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The ratio of rms error 5 to aperture diameter D can be obtained as a function 
of the ground resolution as follows: 


Let R be the range (36,000 km) and r be the resolution In kilometers. Then 



Setting A = 1005 and solving for 5/D yields 
^ = 0.214 x 10" 6 r 

For example, for a ground resolution of 20 km, the value of 5/D = 4.3 x 10’ 6 . 
For a 20 m aperture, 5 = 85 micrometers. On the other hand, for a resolution 
of 6 km and an aperture of 10 m, then 5 = 12.8 micrometers. 

Clearly, very high accuracies will be demanded from the supporting structure 
for the high-frequency radlometry missions In Figure 1. 


2.3 SHAPE CONTROL 

Some shape control is likely to be needed to obtain the required surface 
accuracies. Initial trimming In orbit will probably be desirable. If only to 
reduce the expense of testing before flight. Also, provision should be made 
to adjust the antenna figure to cope with long-term changes In the materials 
due to exposure. 

A worthwhile objective will be to make the structure still enough and ther- 
mally stable enough that It can handle all the short-term excitations without 
deforming too much. Then the shape control system can be of the updating type 
and much less expensive than a full authority system would be. 


2.4 INFLUENCE OF BEAM STEERING AND BAND SWITCHING ON STRUCTURAL REQUIREMENTS 

The radiometer must be able to direct Its beam to any part of the earth's 
disk; thus. It needs to scan about 8 degrees off axis. In addition, the scan 
must be rigid; In order to achieve the desired frequency of coverage, a scan- 
ning rate of hundreds of degrees per minute Is needed. This will cause unac- 
ceptable shaking of the spacecraft If the scan Is entirely mechanical. 
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Therefore, the beam steering will need to be achieved mostly by electronic 
scanning. The slmpllst way to accomplish this is to move the effective feed 
point by varying the gain on individual feed elements (horns, perhaps) in a 
multi-element feed array. Unfortunately, steering the beam by feed movement 
results in large errors for angles more than ten to twenty beam diameters off 
axis. In order to cover the earth, nearly 1,000 beam diameters need to be 
scanned. 

Of course, the art of antenna engineering Is able to achieve much smaller 
errors. One approach, for example, is to design subreflector and reflector 
geometries so as to minimize errors during scan. Another approach that shows 
promise is to use a phased array to illuminate the subreflector. Another pos- 
sibility is to scan rapidly electrically in one direction while slowly moving 
the entire antenna mechanically in the perpendicular direction to cover the 
desired area. 

In addition to steering, the several frequency bands also must be examined. 
The frequency range from 30 to 220 GHz Is obviously too much to be handled 
with a single feed system. Multiple systems will be required, and their loca- 
tion will pose severe problems, especially since they will have to be large in 
order to produce the ±8 degree scan. 

Beam and frequency agility Is the responsibl 1 ty of the antenna engineer. From 
the structural point of view, the need for low spacecraft excitation also 
implies that the dynamic loads on the antenna reflectors will be low. It 
might be possible to ease the beam steering problem by actively shaping the 
subreflector and/or the primary reflector. The amounts of displacements 
required to eliminate the path length error are estimated later In the paper. 

Finally, provision of the needed beam steering with multiple feed systems may 
result In new geometrical configurations for which new structural concepts 
will be required. 

2.5 PACKAGING 

The microwave radiometer operates In geosynchronous orbit. For the purposes 
of this paper, the assumption is made that the deployment will be In geosyn- 
chronous orbit and therefore remote. The launch system Is assumed to be 
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either the Titan IV or a Shuttle-OTV combination, with cargo-bay diameter of 
4.5 meters and an available length of over 10 meters. The Delta launch 
vehicle, with Its smaller launch volume and lower payload, appears to be inap- 
plicable for the ESGP mission. 
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SECTION 3 


DEPLOYABLE STRUCTURAL CONCEPTS 

As in the preceding section, It is convenient to discuss concepts for the 
reflecting surface first. 

3.1 REFLECTING SURFACE 

The requirement that the reflecting surface be solid limits consideration to 
reflective membranes and panels. 

Membrane surfaces can be excellent radio-frequency reflectors. Some care must 
be taken to ensure that the surface conductivity of metallzed plastic films is 
not degraded by cracks in the conductive layer caused by frequent creases. 
For dish-type surfaces, the membrane requires a transverse pressure loading to 
create a wrinkle-free surface of the correct shape. No suitable reflector 
membrane material has low enough in-surface stiffness to enable needed changes 
in Gaussian curvature without Incurring high stresses. Membranes are pliable 
and can be stowed compactly. 

The most prevalent approach for providing a reflector surface Is to use an 
assemblage of stiff panels. A variety of shapes have been proposed, ranging 
from near-hexagonal segments, through ring sectors, to petals. In all cases, 
the panels butt together to produce the large dish-type reflector. For 
launch, the panels are folded or Interleaved to fit In the launch vehicle. 
Each panel Is considered to be stiff and precise enough to maintain Its own 
shape. Panels can be built In several ways, the chief ones being as a honey- 
comb sandwich or a monocoque stiffened shell. 

A novel approach that has been suggested by Composite Optics, Inc. of San 
Diego utilizes a reflector surface composed of a thin flexible shell of 
graphite-epoxy composite. Large areas of the shell can be rolled up for 
launch and allowed to unroll In orbit against a supporting truss structure. 
The shell could comprise the entire surface for smaller antennas. Rolled-up 
shell segments could be stowed with the folded truss for larger apertures. 
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3.2 SUPPORTING STRUCTURE 


Antenna reflector configurations using membrane reflector surfaces must pro- 
vide some means for creating a pressure-type loading across the surface. If 
gas pressure Is used, the reflecting surface is usually joined to a synmetrl- 
cally shaped transparent film around the rim to create a closed pressure 
vessel. Examples are shown in Figure 3 taken from Reference 3. The rim must 
be capable of carrying the compression loading caused by the membranes' pull- 
ing Inward at the rim. The rim may consist of an inflatable torus. The 
assembly is deployed by slow Inflation and is Intuitively very reliable. 

Leakage caused by meteoroid penetration would necessitate a large supply of 
make-up pressurant for long time operation. This can be avoided by making the 
membrane stiff enough to provide its own structural integrity after deploy- 
ment. The ECHO passive satellite, a 100-foot-diameter balloon, was launched 
early in the space age. Its shell was composed of a thin sandwich with Mylar- 
film face sheets and an aluminum-foil core. More recently, technology work In 
Europe has been underway since the early 1980s developing a Kevlar-epoxy com- 
posite surface which Is cured and hardened on orbit after Inflation. See 
References 3 and 4. 

Inflatable antennas, while being vigorously promoted for the lower frequen- 
cies, are generally viewed as being Inapplicable for the high frequencies 
being considered herein. Even when extreme care Is exercised during fabrica- 
tion, the available suitable materials lack the long-term dimensional stabil- 
ity and super-low coefficient of thermal expansion needed for very high 
precision. In addition, Inflatable antennas, once fabricated, are difficult 
to “tune up," even during ground testing. Adjustments In orbit seem to be 
Impossible. 

Membrane antennas shaped and adjusted by electrostatic forces have been pro- 
posed and studied during the last decade. This technique shows good promise 
of being useful, particularly for shallow dishes. Deep dishes are less amena- 
ble to this approach because the high In-surface stiffness of the doubly 
curved membrane causes the shaping pressures to be large. Even for shallow 
dishes, the necessary electrostatic drivers and their support structure tend 
to be heavy and the charged devices must be shielded against arcing due to the 
In-space plasma. On the other hand, rapid adjustment of the lightweight film 
reflector can be accomplished with little disturbance of the spacecraft. 
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Supporting structures for panel -type reflector surfaces are often integrated 
with the reflecting surface itself. Indeed, this approach Is used for the 
many solid dishes flying on communication satellites. Its simplicity is 
attractive and the resulting structure can be made dimensionally stable enough 
to be used at extremely high frequencies. Ingenious concepts have been 
devised for deploying large dishes by hinging between adjacent segments. One 
such technique, termed Sunflower, consists of petals which fold up around the 
symmetry axis and form a complete dish on deployment. The concept underwent 
significant development in the 1960s as a space-borne solar-energy collector 
(see Reference 5). A recently designed descendant Is shown In Figure 4. Note 
that this version deploys to a 15-meter diameter. 

Another segmented-panel approach with Integrated structure was designed for 
high-frequency antennas and Is discussed In Reference 6. Figure 5, taken from 
Reference 6, shows the stack of stowed hexagonal panels, each one of which Is 
rotated Into position and fastened to Its neighbor. Not shown are the mecha- 
nisms required to deploy and attach the segments together. 

Integrated-structure, or panel -only, concepts are attractive because of their 
relative simplicity. They also use well established fabrication techniques 
and appear to be of low risk. They are, however, structurally "thin," so that 
small errors in individual parts grow Into large distortions for large sizes. 
In addition, such structures are difficult to test in a one-g environment. 
Their flexibility combines with the gravity loading to produce deflections 
that are large In comparison to those acceptable for the present application. 
It Is therefore difficult to achieve the desired accuracy, either by fabricat- 
ing the component parts with enough precision or by "trimming" the structure 
by adjustments based on measurements obtained during ground testing. 

The experience and Information obtained by studies and tests over the past two 
decades have shown that structural configurations that are "deep" are much 
more suitable for large high-precision surfaces than are the "thin" ones. See 
References 1, 7, 8, 9, and 10. Not only Is this notion intuitively obvious, 
but also detailed analyses have shown that very high precision is achieved 
with careful fabrication. For example, a recent simulation of a 20-meter- 
diameter tetrahedral -truss structure constructed from 2-meter struts which 
have random lengths with an rms variation of 20 micrometers, showed an 
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Figure 5. View showing the first panel rotated out with its tips displayed and 
lowered for lockup on the center hexagonal hub. 
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expected rms surface error of 43 micrometers. The worst of 100 cases had an 
rms surface error of 72 micrometers. Furthermore, analysis of the deflections 
caused by testing in a one-g field showed an rms error of about 100 micro- 
meters; gravity compensation should be able to decrease that by an order of 
magnitude. 

One antenna with a deep-truss support structure flew on the SEASAT spacecraft. 
As shown In Figure 6, the synthetic aperture radar antenna, which is 10.75 
meters long. Is supported by a deployable truss. The radiating panels are 
stowed and deployed with the truss as seen in Figure 7. This structure, which 
supported an L-band antenna (A * 20 cm) was accurate to better than 2.5 mm 
maximum deflection. This was achieved, and demonstrated with care but without 
heroic efforts; the robustness of the configuration simplified analysis, inte- 
gration and testing. Similar deployment truss concepts have been studied for 
possible use with dish-type reflector antennas, one of these is shown in 
Figure 8. This arrangement has the advantage that It allows the panel segment 
to nest, thereby saving package volume. 

The structural perfoonance of a petal-type deployable reflector can be greatly 
Improved by mounting each petal on a stiffening truss. The approach has been 
suggested by Dornler and is shown in Figure 9. The application is an 8-meter 
reflector for infra-red astronomy. Also being studied for this mission is a 
segmented three-section mirror in which the outer two segments fold inward 
over the central one to form an 8-meter-long package with a 4-meter cross 
section. 

The foregoing truss-stiffened concepts are useful only for diameters smaller 
than the available package length, say up to 10 to 15 meters, depending on the 
launch vehicle configuration. For larger dishes. It will be necessary to 
divide the reflector surface In both directions In the surface. This poses a 
severe problem because almost certainly the surface will have to be cut Into 
segments and stowed separately. The supporting truss can be stowed separately 
also, and the panels can be assembled to the deployed truss by a robot as 
shown In Figure 10. Research is In progress at Langley Research Center on 
such robotic assembly. One concept for the deployable truss which Is being 
extensively studied for various high-precision applications Is the Pactruss 
shown In Figure 11. The deploying truss In this concept Is very strongly syn- 
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Figure 8a. Synchronously deployable Concept B (CREST) 
for stiff-panel reflectors. 
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Figure 8b. Synchronously deployable Concept B (CREST) 
for stiff-panel reflectors. 
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Figure 9. Deployable reflector for FIRST (Dornier System). 
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Figure 1C. Automated curved surface construction concept. 
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chronized and offer's reliable deployment with a few actuators. See Refer- 
ence 11 for the description of recent evaluations of precision application. 


Another concept for constructing large segmented reflectors in remote loca- 
tions is shown in Figure 12. Here, individual modules, each consisting of a 
panel and its associated support truss section (see Figure 13), are stowed in 
a deployment canister which walks around the dish, deploying modules and lock- 
ing each to its neighbors. The development of this Intelligent canister would 
require some effort but seems to be easier than using a robot. Use would be 
made of the fact that each module would be hinged, so far as possible, to its 
neighbors. The hinging would aid In control of the canister motions. 

The furlable, thin-shell reflector panel described in a foregoing section 
might be stowable along with the deployable truss. The rolled-up segments 
could possibly be released after truss deployment and would then settle Into 
frames created by the truss. In this case, the square form of Pactruss would 
probably be more attractive. The panels would then be nearly square. See 
Figure 14. 
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Figure 1 3. The basic truss module. 
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Elevation View 


Perspective View 


Figure 1 4 . Views of Pactruss for offset paraboloid. 




tro Aerospace Corporation • 85-842 


AAC-TN1154 • 1-12-89 


26 


SECTION 4 


PANEL DESIGN 

In the fabrication of panels for large precise antenna reflectors, a mandrel 
Is needed for laying up the panels. The mandrel would either be as large as 
the radius of the paraboloid, or made In several pieces. In either case, the 
expense of the mandrel will be large. 

One way to reduce cost is to make only a few mandrels (one, If possible) and 
replicate panels off of each, using them In the best way to minimize the shape 
error. The following analysis Is aimed at finding the best single mandrel 
shape to produce Identical panels which, when mounted on the support truss at 
the optimum orientation and position relative to the exact paraboloidal 
surface. 

Consider a paraboloid with focal length F, with Its axis along the z axis and 
Its vertex at the origin. Its equation Is 


where 

Jl 2 

r * v x + y 

Let a be the offset of the center of the aperture from the axis of the 
paraboloid and D be the diameter of the aperture. Let p and w be polar coor- 
dinates based on the center of the aperture so that Inside the aperture, 
p < D/2, and « Is measured from the direction of the offset. Then 


r = / a^ + p^ + 2ap cos w 

Consider a circular panel whose center Is located at the location (r 0 , z 0 ) on 
the paraboloid. Let £,n»t be a rlght-and coordinate system, with £ and n tan- 
gent to the paraboloidal surface and C normal to It. Let £ point In the 
meridional direction at the panel center. 
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Let the xz plane pass through the center of the panel. Then 
x = r 0 + £ cos * 0 - C sin «J> 0 
y * n 


z = z 0 + £ sin (ji 0 + cos 


where 


tan 0 O « ~ 


r o 

* If 


Substituting x, y, and z Into the equation of the paraboloid and solving for C 
gives 


C - C* 


where 


£ 2 cos $ 0 + n 2 


+ £ p sin * 0 cos * 0 + 


cos 2 * 0 


+ 4F£ sin 00 - n 2 Sin 2 $ 0 


Let the panel have curvatures in the meridional and circumferential directions 
of k m and k c , respectively. Also let the center displacement of the panel In 
the C direction be Co and the tilt In the meridional plane be a. The equation 
of the panel surface Is then 


t * tp + Co + a£ 


where 


Cp * j( k mC^ + 


Then the local error In the normal direction between the panel and the para- 
boloidal surface Is 

& » C p - C* + Co + a£ 
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The mean-squared error Is given by 


S 


2 

rms 


y c * - c p - 



- 



where the integrations are carried out over the area of the panel. The mean- 
squared error is minimized when 


Co 


a 



d£dn 



£d£dn 


This process yields the rms error for a particular value of r. The mean- 
square error for the entire antenna Is obtained by averaging over the 
aperture. 

The computer code UNIPANL.C, which Is Included In Appendix A, was written to 
perform the indicated integrations and averages, and determine the rms error 
for the antenna. The program is Interactive, requesting Inputs fo D, F, and 
offset, then repeatedly asking for the panel size and ratio of circumferential 
to meridional curvature. The Integrations are performed numerically with five 
intervals in the radius and 15 degree Intervals around the circumference. The 
panel curvature that gives the least rms error over the entire aperture is 
found by a stepping type of search for the minimum. 

Some results for panels which have the same curvature In both directions 
(spherical mandrel) are shown In Figure 15. Note that using an offset feed 
with an F/D of 1.5 yields almost the same results for Inaccuracy as those for 
a centered-feed antenna with F/D * 1.0. To understand these results, consider 
a 20-meter diameter to be used at a frequency of 100 GHz and require A/100 
accuracy. Then Sp^/D * 1.5 x 10"®. With a centered-feed and F/D * 1.5, the 
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Figure 15. Antenna surface error caused by identical spherical panels. 
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panel size could be as large as 2 meters. For an offset feed and F/D of 1.5, 
the allowable panel size Is only 1 meter. Note also that a resolution of 6 km 
for the ESGP radiometer would need a value of S/D of about 1.3 x 10" 6 . For 
the offset feed case, there would be about 20 panels needed across the aper- 
ture diameter. 

Incidentally, some trials with the circumferential curvature slightly higher 
than the meridional Indicates significant reduction in the error. Also, pro- 
viding two mandrels would help a great deal. 
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SECTION 5 


VARIABLE GEOMETRY 

One approach to avoiding pattern deterioration when scanning would be to 
adjust the shape of the reflector as the scanning occurs. In order to deter- 
mine the magnitude of the motions required of the surface, an analysis was 
made of the path-length error due to scanning. The approach used was to find 
the tilted paraboloid for which the mean-square normal distance from the orig- 
inal paraboloid was a minimum. This analysis Is coded In the program ADJUST. C 
Included In the appendix. 

A sample of the output of ADJUST Is Included In Table 1. The case treated is 
a 20-meter-dlameter offset-feed antenna with an F/D * 1.5 and an offset of 
12.5 meters. The rms value of the correction is about 1.5 cm and the maximum 
value Is about 5 cm. These are sizable motions, but not nearly as large as 
would occur If the beam were steered by rotating the entire antenna. 

The indicated surface adjustment would be accomplished by actuators. If the 
surface were a continuous one, say an electrostatically controlled membrane, 
then the surface would tend to fair the shape between control points. If the 
surface Is made up of segmented panels, then the control would be applied at 
the attachment points. Since the panels would each be shaped to conform to 
the untilted paraboloid, they would exhibit some unavoidable residual error 
when trying to fit the scanned paraboloid. The program ADJUST Includes the 
ability to examine Individual panels for their residual errors. Results for 
the worst-case panels are shown In Figure 16. Examination shows that the 
residual errors are similar to those due to using Identical panels. 

The foregoing results are calculated for scanning by simple feed motion. Much 
smaller errors will result from the more advanced scanning techniques that 
will be used. If variable geometry Is used, the motions and residual errors 
would be accordingly smaller. 
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TABLE 1: REFLECTOR CORRECTIONS FOR SCAN 


F = 30.000000 D = 20.000000 rO = 12.500000 delta = 8.000000 
Displacement of focal point = -3.916113, 0.000000, -1.946797 
Rms path length error = 0.027682 

New focal length = 27.796289 Rms correction = 0.015121 



0.2 

0.4 

0.6 

0.8 

0 

-1 . 481e-003 

-6. 442e-003 

-1 . 562e-002 

-2. 967e-002 

15 

-1 . 307e-003 

-5. 733e-003 

-1.400e-002 

-2.677e-002 

30 

-8. 270e-004 

-3.783e-003 

-9. 546e-003 

-1.875e-002 

45 

-1.670e-004 

-1 . 080e-003 

-3. 332e-003 

-7.486e-003 

60 

5 . 027e-004 

1 . 701e-003 

3. 140e-003 

4. 380e-003 

75 

1 . 009e-003 

3.867e-003 

8. 314e-003 

1.409e-002 

90 

1 . 223e-003 

4. 888e-003 

1 . 098e-002 

1 . 946e— 002 

105 

1 . 093e-003 

4.538e-003 

1 . 057e-002 

1.943e-002 

120 

6. 600e-004 

2.957e-003 

7 . 369e-003 

1.437e-002 

135 

4 . 544e-005 

6. 180e-004 

2 . 386e-003 

6.028e-003 

150 

-5.790e-004 

-1 . 801e-003 

-2. 867e-003 

-2. 950e-003 

165 

-1 . 040e-003 

-3.602e-003 

-6. 814e-003 

-9.766e-003 

180 

-1 . 209e-003 

-4 . 266e-003 

-8.276e-003 

-1.230e-002 

195 

-1 . 040e-003 

-3. 602e-003 

-6. 814e-003 

-9.766e-003 

210 

-5 . 790e-004 

-1 . 801e-003 

-2.867e-003 

-2. 950e-003 

225 

4 . 544e-005 

6. 180e-004 

2 . 386e-003 

6. 028e-003 

240 

6.600e-004 

2. 957e-003 

7 . 369e-003 

1.437e-002 

255 

1 . 093e-003 

4. 538e-003 

1 . 057e-002 

1. 943e-002 

270 

1 . 223e-003 

4 . 888e-003 

1 . 098e-002 

1 . 946e-002 

285 

1 . 009e-C03 

3.867e-003 

8.314e-003 

1.409e-002 

300 

5 . 027e-004 

1. 701e-003 

3 . 140e-003 

4. 380e-003 

315 

-1 . 670e-004 

-1.080e-003 

-3.332e-003 

-7.486e-003 

330 

-8. 270e-004 

-3. 783e-003 

-9.546e-003 

-1. 875e-002 

345 

-1.307e-003 

-5.733e-003 

-1.400e-002 

-2. 677e-002 


psi = 0.000000 


1.0 

-4 . 924e-002 
-4.466e-002 
-3. 198e-002 
-1.407e-002 
5 . 004e-003 
2. 095e-002 
3. 030e-002 
3. 131e-002 
2.442e-002 
1.223e-002 
-1 . 202e-003 
-1. 151e-002 
-1 . 536e-002 
-1 . 151e-002 
-1 . 202e-003 
1.223e-002 - 
2.442e-002 
3. 131e-002 
3 . 030e-002 
2.095e-002 
5.004e-003 
-1 . 407e-002 
-3 . 198e-002 
-4 . 466e-002 
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SECTION 6 


CONCLUDING REMARKS 

The study reported herein Is only a beginning. There remains a great deal of 

Investigation before a good configuration can be selected for development. 

Among the questions are: 

• What rms accuracy Is needed by the radlometry mission? A/30? A/50? 

A/100? 

• How good can electronic scanning be? Feed-motion scanning becomes unac- 
ceptable at 10 beamwldths. The mission needs 1,000. 

• Can robots or Intelligent canisters be developed In time to be available 
for remote assembly of antennas needed In the year 2000? 

• Can long-time microstrain stability for the available materials be 
assured? 

• What are the magnitude and distribution of the forces required to adjust 
the shape of continuous reflecting surfaces? 

• How accurately can large continuous shells be built? 

The future Is promising. 
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PRECEDING PAGE BLANK NOT FILMED 



*** SCITECH. H *** 


Date of file: 12/ 5/88 


Printed on: 12/05/88 


* SCITECH. H - Header file for all science, engineering, and technical 

* programs. This "included all the needed headers. 

* J. M. Hedgepeth 7/31/86 


* 


* 

Added 

string. h 

and process. h to list. 

jmh 

12/30/86 

♦ 

Added 

float . h. 


jmh 

2/10/87 

* 

Added 

io. h 


jmh 

9/4/87 

* 

* 

*/ 

Added 

macro for 

flushconO and constants 

jmh 

11/23/88 


#include 

♦include 

♦include 

#include 

♦include 

♦include 

♦include 

♦include 

♦include 

♦include 

♦include 

♦include 

♦include 

♦include 

♦include 


' stdio. h> 
<conio. h^ 
<io.h> 
<ctype. h> 
<dos.h> 
<malloc. h> 
<math.h> 
<float.h> 
<errno. h> 
<stdlib . h> 
<string.h> 
<t ime. h> 
<scieng.h> 
<fcnt 1 . h> 
<process.h> 


/* Force use of functions, not macros 
*/ 

♦undef toupper 
♦undef tolower 


/*********** 

*/ 

♦define flushcon() while (kbhit ( ) ) getch() 

/* Values of particular mathematical constants 
*/ 

♦define PI 3.14159265358979 

♦define EULER 0 . 577215664901532 

♦define EBASE 2.718281828459045 

/**************«♦***♦***♦♦*♦«*♦♦**♦♦♦♦♦*♦♦♦**♦**♦♦ 

*/ 


40 



*** UNIPANL.C *** 


Date of file: 11/10/88 


Printed on: 01 09 89 


/* UNIPANL.C - Determine errors in surface shape if a paraboloid is approximate 

* by panels of uniform meridional and circumferential curvature. 

* 

* J. M. Hedgepeth 11/9/88 

* 

*/ 

#include scitech.h) 

#define flushconH whilefkbhit ( ) ) getch() 

double F, rO,D, curv, ratio,d panl; 
char line [80] ; 

/* 

* F 

* rO 

* D 

* curv 

* ratio 

* d_panl 

* 

*/ 

double degtorad; 

double f indmin(double [ ] ) , getsumsq(double, double []), panel_err (double, double 

mainf ) 

{ 

int i,j,k; 

double xv , yv , zv , xf , yf , zf , tempO , tempi , temp2 , rms , path, *ptr ; 
static double eps[120]; 
char chr; 

degtorad = (double)PI/180. ; 
setmode( fileno(stdprn) ,0_TEXT) ; 

while(l) { 

printf ( tf \n\nEnter the reflector focal length, aperture diameter , \n" 

"and distance to aperture center. \n? "); 
gets(line) ; 

for(i=0; i<80 8A line[i]; i++) 
if (line[i] == *,*) 
line[ i] = * * ; 

if (sscanf ( line, tf *lf %lf %lf", &F, &D, &r0) == 3) 
break; 

printf ("\n\007Incorrect entry. " ) ; 

\ 

j 

while(l) { 

while(l) { 

printf( M \n\nEnter panel diameter and ratio of circumferential to 


Focal length of base paraboloid 

Distance from paraboloid axis to aperture center 

Aperture diameter 

Meridional panel curvature 

Ratio of panel circumferential to meridional curvature 
Panel diameter 
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"mer idianal curvatures. n° 


*t . 


gets ( 1 ine ; 

for(.i-0; i<80 && linefij; 
if(line[i] == 1 , ’ ) 
line[i] = ’ ’ ; 

if (sscanf( line, ” Slf Slf", &d_panl, &ratio) == 2) 
break; 

printf ( "\n\007Incorrect entry. " ) : 


nns - f indmin ( eps ) ; 

printf ("\n\n REFLECTOR ACCURACY FOR PANELS WITH UNIFORM " 

"CURV ATURE\n\n" ) ; 

printf ( "F = Slf D - Slf rO = Slf d = ^lf ratio = *lf\n", 

F, D, rO.dpanl, ratio) ; 

printf ("Meridional curvature for minimum rms surface error = °ilf\n", 

curv) ; 

printf ("Rms surface error of reflector = *14. 31e\n\n" , rms) ; 
printf ( "Press any key to continue? "); 
f lushcon( ) ; 
getch( ) ; 

printf("\n 0.1 0.3 0.5 " 

"0.7 0.9"); 

for(i=0; i < 1 3 ; i++) { 

printf ("\n %3d ",i*15); 

for ( j=0 ; j<5; j++) 

printf("Ssl4.31e",eps[5*i + j ] ) ; 

} 

printf (”\nDo you want to print the results? <N> "); 

if ( toupper(getche( ) ) == ’ Y* ) { 

fprintf (stdprn, "\n\n\n\n REFLECTOR ACCURACY FOR " 

"PANELS WITH UNIFORM CURVATURE\n\n") ; 
fprintf (stdprn, "F = %lf D = *lf rO = %lf d = Slf ratio = Slf" 

"\n",F,D,r0,d_panl, ratio) ; 

fprintf (stdprn, "Meridional curvature for minimum rms surface error 

”= Slf\n",curv) ; 

fprintf (stdprn, "Rms surface error of reflector = S14.31e\n\n" , rms) 
fprintf (stdprn, " 0.1 0.3 0.5" 

" 0.7 0.9"); 

for(i=0; i<13; i++) { 

fprintf (stdprn, "\n S3d ",i*15); 
for( j=0; j <5 ; j++) 

fprintf ( stdprn, "S14.31e",eps[5*i + j]); 

} 

fprintf (stdprn, "\n\f") ; 
fflush( stdprn) ; 
fflush(stdin) ; 


printf ("\n\nDo you want to examine new panel parameters? <N> "); 

f lushcon( ) ; 

if(toupper(getche( ) ) != *Y* ) 
break; 

} 

putchar( ’ \n’ ) ; 
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void flushcon' 

whi le( kbhit ; ; i 
getch( ) ; 

i 


/t*******t********************************************* 

%/ 


43 



Izsq - Iz ~ Izksi = Iksisq - I - 0; 
phi - at.an< r 2. . F ) ; 
cosphi = cos (phi}; 
sinphi = sin(phi); 

for(i=0; i < 12 ; i++) { 

omega - { (double) i + 0.5)*15. *degtorad; 
cosom = cos (omega); 
sinom = sin(omega); 
for j -0 : j'5; j++) { 

rho = d_panl*(2*j + l)/20.; 
ksi = rho*cosom; 
eta - rho*sinom; 

B - 2»*F/cosphi + ksi*sinphi*cosphi ; 

C - ksi*ksi*cosphi*cosphi + eta*eta; 
zeta = C/(B + sqrt(B*B - C*sinphi*sinphi) ) ; 
zeta - = kx*(ksi*ksi + ratio*eta*eta)/2. ; 
eps [ i ] [ j] = zeta; 

I += rho; 

Iz += rho*zeta; 

Izsq += rho*zeta*zeta; 

Izksi + = rho*zeta*ksi ; 

Iksisq += rho*ksi*ksi; 

/ 

} 

dave = Iz/I; 
sx = Izksi/I; 

rms = sqrt(Izsq/I - davetdave -sx*sx*Iksisq/I) ; 
if(kbhit()) f 

printf("\n\n PANEL ERRORS\n\n") ; 

printf ("F = *lf D = *lf rO = %lf d = *lf ratio = \lf\n\n", 

F, D, rO.dpanl, ratio) ; 

printf ( "Meridional curvature = *lf Panel location = *lf\n\n", kx, r) 
printf ("Rms surface error of panel = %14.31e\n",rms) ; 
printf("\n 0.1 0.3 0.5 

"0.7 0.9\n”) ; 

for(i=0; i<12; i++) { 

omega = ( (double) i + 0.5)*15.; 
printf ( "\n*6. Ilf" , omega) ; 
omega *= degtorad; 
cosom = cos (omega); 
for(j=0; j<5; j++) { 

rho = d_panl*(2*j + l)/20.; 
ksi = rhoftcosom; 

zeta = eps[i][j] - dave -sxtksi; 
printf ("\14. 31e" , zeta) ; 

} 

} 

printf ("\n\nPress any key to continue..."); 
flushcon( ) ; 
getch( ) ; 

} 

return rms; 


/**************************************************************** 
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************************************************************ 

*/ 

double findmint double eps[]) 
int n; 

double temp, interval, sumsq; 
curv = 0.5/F; 

sumsq = getsumsq(curv,eps) ; 
interval = -0.01/F; 

for'n=0; n< 10; n++) { 

print f f M \nCycle %2d Interval = *lf\nCurvature, rms = %lf, °slf", 

n, interval , curv, sqrt ( sumsq) ; 
while((temp = getsumsq( curv + interval , eps ) ) <= sumsq) { 
curv += interval; 
sumsq = temp; 

printf (" *lf, Slf", curv, sqrt ( temp) ) ; 

\ 

j 

interval /= 2.; 

i 

/ 

return sqrt ( sumsq) ; 

i 

/t******************************************************** 

*/ 

double get sumsq (double cur, double eps[]) 

{ 

int i,j,k; 

double temp, sum, omeg, cosom, r , factor; 

static double weight[] = { 0.1, 0.3, 0.5, 0.7, 0.9}; 

for( i=0, sum=0; i < 1 3 ; i++) { 
factor = 1. ; 
if ( i == 0 ! i == 12) 
factor :: 0.5; 
omeg = i * 15 . *degtorad; 
cosom = cos (omeg); 
for(j=0; j<5; j++) { 

temp = (2*j + l)*D/20.; 

r = sqrt(r0*r0 + 2 . *temp*r0*cosom + temp*temp); 
k = 5*i + j; 

eps[k] = panel_err(cur,r) ; 

sum +— factor*weight[ j]*eps[k]*eps[k] ; 

} 

} 

return sum/30; 


/^tt****t***********t******t********************** 

*/ 

double panel_err(double kx, double r) 

{ 

int i , j , k, lp[3] ; 

double B , C , ksi , eta, zeta, rho, omega, dave, sx, rms , Izsq, Iz , Izksi , I , Iksisq, phi ; 
double cosom, sinom, cosphi , sinphi ; 
static double eps [12] [5]; 
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*** ADJUST.C *** 


Date of file: 10/ 9/88 


Printed on: 01/09/89 


/* ADJUST. 
* 

* 

* 

* 

* 

* 

* 

* 

*/ 


- Determine correct paraboloid for scanned beam and correction 
distances normal to unscanned paraboloid. 

J. M. Hedgepeth 9/19/88 

Added minimization of sum-displacements-squared with 

respect to F. 9/27/88 

Added examination of residual errors for panels. 9/29/88 

Changed to best-fit panel motion. 9/30/88 


♦include scitech. h> 

double F, Fnew, delta, psi , rO , D , a,b , c; 


/* 

* 

* 

♦ 

♦ 

♦ 

♦ 

♦ 

* 

*/ 

double cosdel , sindel , cospsi , sinpsi , degtorad; 

double normdist (double, double) , findmin( double []) ,getsumsq( double, double []); 

double minpath(double, double, double); 

void panel_err(double []); 

raain( ) 

{ 

int i,j,k; 

double xv, yv, zv,xf ,yf , zf , tempO, tempi , temp2, rms , path, *ptr; 

static double eps[120]; 

char chr; 

degtorad = (double)PI/180. ; 

setmode(f ileno(stdprn) ,0_TEXT) ; 

while ( 1) { 

printf ( M \n\nEnter the reflector focal length, aperture diameter, \n" 
"and distance to aperture center. \n? "); 

gets( line) ; 

for(i=0; i<80 && 1 ine [ i ] ; i++) 
if ( line[i] == V) 
linefi] = ’ 

if (sscanf ( line, " *lf *lf *lf", &F, &D, &r0) == 3) 
break; 


r line [80] ; 


F Focal length of base paraboloid 

Fnew Focal length of scanned paraboloid 

delta Angle of scan in degrees 
psi Azimuth of scan in degrees 

rO Distance from paraboloid axis to aperture center 

D Aperture diameter 

a, b, c Parameters of scanned paraboloid 
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print f ” n 007 Incorrect entry. 


while(l) { 

whiled) { 

printf ("\n\nEnter scan angle and azimuth in degrees. \n? 
gets ( line) ; 

for(i=0; i<80 && line[i]; i++) 
if( 1 ine[ i] == ’ , * ) 
line[i] = * * ; 

if i. sscarif ( line, " Xlf Xlf", &delta, &psi) == 2) 
break ; 

printf (dnXOOTIncorrect entry. ") ; 

i 

s 

cosdel = cos(degtorad*delta) ; 
sindel = sin(degtorad*delta) ; 
cospsi = cos(degtorad*psi) ; 
sinpsi = sin(degtorad*psi) ; 

rms = findm in(eps) ; 

xv = rO - a*cospsi*cosdel - b*sinpsi - c*cospsi*sindel/Fnew; 
yv = -a*sinpsi*cosdel + b*cospsi - c*sinpsi*sindel/Fnew; 
zv = rO*rO/4./F + a*sindel - c*cosdel/Fnew; 

xf = xv + Fnew*cospsi*sindel; 
y f = yv + Fnew*sinpsi*sindel; 
zf = zv + Fnew*cosdel - F; 

path = minpath(xf, yf, zf + F); 

printf ("\n\n REFLECTOR CORRECTIONS FOR SCAN\n\n"); 

printf ("F = Xlf D = Xlf rO = Xlf delta = Xlf psi = Xlf\n", 

F,D, rO, delta , psi) ; 

printf ( "Displacement of focal point = Xlf, Xlf, Xlf\n ,xf,yf,zf); 
printf ("Rms path length error = Xlf\n",path) ; 

printf ( "New focal length = Xlf Rms correction = Xlf\n\n ,Fnew,rms), 

printf ("Press any key to continue? "); 

getch(); 

printf("\n 0.2 0.4 0.6 

"0.8 1.0"); 

for(i=0; i< 24 ; i++) { 

printf("\n X3d ",i*15); 
for(j=0; j<5; j++) 

printf("X14.31e",eps[5*i + j ] ) ; 

printf ("\nDo you want to print the results? <N> "); 

if ( toupper(getche( ) ) == * Y* ) { 

fprintf (stdprn, "\n\n\n\n REFLECTOR CORRECTIONS FOR 

" SCAN\n\n") ; 

fprintf (stdprn, "F = Xlf D = Xlf rO = Xlf delta = Xlf psi = Xlf" 

"\n", F,D,r0, delta, psi ) ; 

fprintf (stdprn, "Displacement of focal point = Xlf, Xlf, Xlf\n", 

xf,yf,zf); 

fprintf ( stdprn, "Rms path length error = Xlf\n",path) ; 

fprintf (stdprn, "New focal length = Xlf Rms correction = Xlf\n\n", 

Fnew, rms) ; 
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fprintf stdprn," 


0.2 


0.0 

1 . 0 ": ; 


0. I 

0.8 

for ! i “0 : i 24 : 5 ( 

fprintf ( stdprn, n °o3d " , 1*151 ; 
for(j=0; j < 5 ; j + + ' 

fprintf stdprn, " , il4.31e",epsf5*i j]); 

} 

fprintf (stdprn, "\.n' f" ' : 
fflush( stdprn) ; 

t 


printf ( "\n nDo you want to examine individual panels? <N> ") 

if ( toupper 'getchei ) ' == * Y’ ) 

panel errieps ) ; 

while(i) { 

printf ( "\n\nEnter ’ F* to change focal length, (%10.31f)" 
"\n ’D* to change diameter, (%10.31f)" 

"\n ’0’ to change offset, (%10.31f)" 

”\n ’q’ to quit," 

"\n <:CR> to continue." 

"\n? " , F, D, rO) ; 

if((chr = toupper(getche( ) ) ) == ’Q’) 
ex i t ( 0 ) ; 

if (chr == ’F’ I! chr ==’D’ :: chr == ’O’) { 
printf ("\n\nEnter new value? "); 
gets( line) ; 

if (sscanf ( line, " ^lf" ,&temp2) == 1) { 
if (chr == ’F’) 

F = temp2; 
if (chr == ’D’) 

D = temp2; 
if (chr == ’O’) 
rO = temp2; 


} 


} 

else if(chr == *\r’) 
i = 0; 

else 

putchar( ’ \007* ) ; 


/************************************************************* 

*/ 


double normdist( double x, double y) 

{ 

double r, theta, dif , cosdif .sindif ,cosphi,sinphi; 

double alfal,alfa2,alfa3,betal,beta2,beta3, A, B,C, tempO, tempi; 

r = sqrt(x*x + y*y) ; 
theta = atan2(y,x); 


dif = psi*degtorad - theta; 
sindif = sin(dif ) ; 
cosdif - cos(dif); 


cosphi = cos(atan(r/2./F) ) ; 
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sinphi : cosphi*r 2. , F; 


alfal - cosdif*cosdel*sinphi +■ sindel*cosphi ; 
alfa2 = sindiftsinphi; 

alfa3 = cosdif*sindel*sinphi - cosdeltcosphi ; 

betal = a/F + (cosdif*r - cospsifcrO 1 ' *oosdel/F 

- { r*r - rO*rO)*sindel/4./F T; 
betal! - b/F + (sindif*r - sinpsi*rO)/F; 
beta3 = c/F/Fnev. + (cosdif*r - cospsi*rO)*sindel/F 

+ ( r*r - r0*r0 ) *cosdel/4 . /F/F : 


A = alfal*alfal + alfa2*alfa2; 

R = -F*'alfal*betal + alfa2*beta2) r 2. *Fnew*alfa3; 

C = F*F*(betal*betal + beta2*beta2) - 4. *F*Fnew*beta3; 

terapO = B*B - A*C; 
iff tempO < 0) 

return l.0e99: 
tempO = sqrt(tempO); 

if (B > O' 

return ((-B + tempO)/A); 

else 

return ((-B • tempO)/A); 


/t*****************************^********************* 

*/ 

double f indminfdoubie eps[]) 

r 

y 

Lnt n; 

double ratio, temp, interval , sumsq; 
ratio = 1. ; 

sumsq = getsumsqfratio.eps) ; 
interval = 0.01; 

for(n=0; n<10; n++) { 

printf ( M \nCycle *d ",n); 

whilefftemp = getsumsqf ratio + interval.eps)) <= sumsq) { 
ratio += interval; 
sumsq = temp; 

} 

interval /= -2.; 

} 

Fnew = ratio*F; 
return sqrtf sumsq); 

) 


/ ********************************************************* 
*/ 


double getsumsqf double ratio, double eps[]) 

{ 

int i,j,k; 

double tempO, tempi , sum, omeg, cosom, sinom,x,y, r , theta; 
static double weight[] = { 1.0, 2.0, 3.0, 4.0, 2.5}; 
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Fnew * ratio*F; 


tempO = rO/2./F; 

a = 2. *Fnew*(sindel ■+■ tempO*cospsi*cosdel)/(cosdel - terapO*cospsi*sindei 
b = 2 . *Fnew*tempO*s inps i / ( cosde 1 - tempO*cospsi*sindel) ; 
c = (a*a + b*b)/4. ; 

for( i=0 , sura=0 ; i < 24 ; i++) ( 
omeg - i*15 . *degtorad; 
eosom - cos (omeg); 
sinom - sin(omeg); 
for ( j=0 ; j<5; j^+) { 

tempi = Cj + 1)*D/10.; 

x = rO + tempi *cosom; 

y = templ*sinom; 

k = 5*i + j; 

eps[k] - normdist (x,y) ; 

sum +■= weight [ j ] *eps [k] *eps [k] ; 

} 

i 

j 

return sum/300.; 


*/ 

void panel err (double eps[]) 

{ 

int i , j , k, lp[3] ; 

double x , y , r , xO , yO , rho , omega, tempO , dave , sx , sy , rms ; 
static double weight[] = { 1.0, 2.0, 3.0, 4.0, 2.5}; 
double Mat [12] ; 

while(l) { 

while(2) { 

printf ("\nEnter coordinates (xO, yO) of panel center... ? "); 
gets( line) ; 

for(i-0; i<80 && line[i]; i++) 
i f ( 1 ine [ i ] == ’,*) 
line[i] = * *; 

if(sscanf(line," *lf *lf", &x0,&y0) == 2) 
break; 

else 

putchar ( ’ \007 * ) ; 

} 

while(2) {- 

printf("\nEnter size (diameter) of panel... ?•"); 
gets (line) ; 

if (sscanf (line, " *lf", &r) == 1) 
break; 

else 

putchar ( *\007* ) ; 

J 

r /= 2. ; 

for(i=0; i < 12 ; i++) 

Mat[i] = 0; 

for(i=0,rms = 0; i < 12 ; i++) { 

omega = (double) i*30. *degtorad; 
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for j-0; j 5; j— . 

rho ‘ r*(j r 15.: 
x = rho*cos ( omega ') ; 
y = rho*sim omega; ; 

k = 5*i + j; 

eps[k] = normdist(x0 x,y0 ■“ y); 

Mat [3] += weight ■ j ] ; 

Mat [ 1 ] += weight [j]*x; 

Mat [2] += weight[j]*y; 

Mat [3] += weight [j]*eps[k] ; 

Mat [ 5 j += weight [j]*x*x; 

Mat [6] += weight [j ] *x*y; 

Mat[7] += weight[ j]*x*eps[k] : 

Mat [10] + = weight [ j]*y*y; 

Mat [11] + = weight [ j]*y*eps[k] ; 

} 

Mat [4] = Mat [ 1] ; 

Mat [8] = Mat [2] ; 

Mat [9] = Mat [6] ; 

dcrout (4, 3 , 1 ,Mat , . 000001 , &temp0 ,&i , lp) ; 

dave = Mat [3] ; 
sx = Mat [7] ; 
sy = Mat [ 11] ; 

for(i=0,rms = 0; i < 12 ; i++) { 

omega = (double) i*30.*degtorad; 
for(j=0; j<5; j++) { 
rho = r*( j + l)/5. ; 
x = rho*cos (omega) ; 
y = rho*s in ( omega) ; 
k = 5*i + j; 

epsjk] -= dave + sx*x + sy*y; 
rms += eps[k]*eps[k] ; 

j 

} 

rms = sqrt(rms/150. ) ; 

printf ( ”\n\n RESIDUAL PANEL ERRORS\n\n" ) ; 

printf ("F = %lf D = *lf rO = Xlf delta = %lf psi = *lf\n", 

F, D,rO, delta, psi) ; 

printf ( "Location of panel center = *lf, \lf ",x0,y0); 
printf ("Panel size = *lf\n",2*rl; 

printf("Best fit panel adjustment (dcenter.sx.sy) = *lf, *lf, %lf\n", 

dave,sx,sy); 

printf ("Rms residual error = %lf\n\n" , rms) ; 

printf ("\n 0.2 0.4 0.6 " 

"0.8 1.0"); 

for(i :; 0; i < 12 ; i++) { 

printf (’’\n *3d ",i*30); 

for(j=0; j <5 ; j++) 

printf ("%14.31e",eps[5*i + j]); 

\ 

printf ("\n\nEnter ’q’ to return to main program, \n" 

" Any other key to examine another panel... "); 

if ( toupper( getche( ) ) == ’ Q’) 
return; 
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* 

double minpath(double xf, double yf, double zf) 

1 

int i,j,k; 

double sum, sumsq, x , y, z , r , rho , omega, theta, lengt.h,weight,denom,Z; 

denom - 0 ; 
sum = 0 ; 
sumsq = 0; 

for(i ~ 0; i<24; i++) ( 
omega = i*15*degtorad; 
forf j -0; j<5 ; j++) ( 

rho = ( j + 1)*D/10. ; 
x - rO + rho*cos ( omega ) ; 
y = rho*s in ( omega) ; 
r - sqrt(x*x + y*y) ; 
z = r*r/(4.*F); 
theta = atan2(y,x); 

x - = xf; 

y -= yf; 

length = -zf - z*cosdel - r*cos( theta - psi*degtorad) *sindel + 
sqrt(x*x + y*y + (z - zf)*(z - zf)); 

weight = j + 1; 
if ( j == 4) 

weight /= 2.; 

denom += weight; 

sum += weight*length; 

sumsq += weight*length#length; 

} 

} 

Z = sum/denom; 
sumsq /= denom; 
sumsq -= Z*Z ; 
return sqrt(sumsq); 


ttttttttttttttttttttttt*************************** 

*/ 
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